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ABSTRACT 

We present a Chandra observation of IRAS 19254-7245, a nearby ULIRG also 
known as the Superantennae. The high spatial resolution of Chandra allows us 
to disentangle for the first time the diffuse starburst emission from the embedded 
Compton-thick AGN. The 2-10 keV spectrum of the AGN emission is fitted by a 
flat power-law (r = 1.3) and a He-like Fe Ka line with EW~1.5 keV, consistent 
with previous observations. The Fe Ka line profile could be resolved as a blend 
of a neutral 6.4 keV line and an ionized 6.7 keV (He-like) or 6.9 keV (H-like) line. 
Variability is detected compared with the previous XMM-Newton and Suzaku 
observations, demonstrating the compact size of the iron line emission. We fit 
the spectrum of the galaxy-scale extended emission excluding the AGN and other 
bright point sources with a soft thermal component with kT ~ 0.8 keV. The 
luminosity of the extended emission is about one order of magnitude lower than 
that of the AGN. The basic physical and structural properties of the extended 
emission are fully consistent with a galactic wind being driven by the starburst 
(no contribution to the feedback by the AGN is required). A candidate ultra- 
luminous X-ray source is detected 8" south of the southern nucleus. The 0.3-10 
keV luminosity of this off-nuclear point source is ~ 6 x 10 40 erg s _1 if the emission 
is isotropic and the source is associated with the Superantennae. 

Subject headings: galaxies: active-galaxies: starburst-X-rays: galaxies-galaxies: 
individual: IRAS 19254-7245 (the Superantennae) 



-2- 



Introduction 



Ultra-Luminous Infrared Galaxies (ULIRGs) emit most of their total energy output in 
the far-infrared (FIR) domain, with a FIR luminosity greater than 1O 12 L . Many ULIRGs 
are a composite of an active galactic nucleus (AGN) and a starburst (SB), with the violent 
star formation activity triggered by a galaxy collision or merger. The emission from the cen- 
tral AGN and the dense star-forming regions is obscured by dust and gas, and the fractional 
contributions to the luminosity by the AGN and SB are not well constrained. Obtaining 
better insights into the energetics of these systems is crucial to understand the interplay 
between supermassive black holes (SMBHs) and the formation and evolution of galaxies, 
which is one of the major unsolved problems in astrophysics. 



X-ray observations of ULIRGs abo ye 3 keV, performed with XMM-Ne wton ( IFranceschini et al 



20031 : iBraito et al.ll2003f ) and Chandra flPtak et allbOOa iTeng et al.ll2005f ) have proved to be 
a fundamental tool to disentangle the contribution of SB and AGN activity and to investi- 
gate the presence of hidden AGNs in these sources. Indeed some of the ULIRGs, classified as 
pure SBs based on optical and IR spectroscopy, show s pectral properties typical of obscured 
AGNs only when obse rved in X-rays (i.e., NGC 6240: llwasawa fc Comastrilll998t Arp 299: 
Delia Ceca et al.l 120021 ) . Previous X-ray surveys of ULIRGs have also shown that thermal 
emission from a hot plasm a with a rather constant temperature kT ~ 0.7 keV is always 
present ( iGrimes et al.ll2005l ). This emission dominates the X-ray spectra below 1 keV and it 
is likely associated with a nuclear or circumnuclear starburst. In particular, Chandra obser- 
vations of luminous and ultraluminous IR galaxies have proved to be essential to investigate 
the spatial structure and physical properties of this soft diffuse X-ray emission. Thanks to 
its high spatial resolution, Chandra observations have shown that in ULIRGs, like in lo- 
cal starbursts, this soft diffuse emission extends far beyond the sta rburst region itself. This 



emission is primary evidence of the presence of ga lactic superwinds ( iChevalier fc Cleggill985 



Heckman et al.lll993l ; IStrickland fc Stevens! 120001 ). which are believed to play a major role in 
the evolution of galaxies and metal enrichment of the intergalactic medium. However, these 
studies have also shown that ULIRGs are intrinsically faint X-ray sources, thus an X-ray 
follow up is feasible only for the very brightest ones. 

Understanding the physical processes at work in ULIRGs is important since observa- 
tional evidence shows that ULIRGs are generally advanced mergers of gas-rich galaxies, which 
are considered to be the origin of many massive elliptical and SO galaxies. If indeed ULIRGs 
are forming spheroids, an AGN/SB association is naturally supported by the evidence that 
all spheroidal galaxies host quasar relics in the form of supermassive black holes (SMBHs), 
and the quasar stage could be a phase during the evolution of these syste ms. Support 



for this view has come from high-resolution spectroscopy of ULIRG mergers (IGenzel et al. 
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200ll ). whose dynamical properties have proved that they indeed are "ellipticals in forma- 



tion" . In this framew ork high resolution simulations of galaxies mergers (IHopkins et al.ll2006 



Springel et al.l 120051 ) have shown that galaxy mergers may play a fundamental role in super- 
massive black hole growth and galaxy evolution. Indeed the merger provides the inflow of gas 
which is likely to fuel both the starburst and, at some stages, an AGN. These simulations 
predict that for most of this phase the AGN is obscured, with the optically-visible AGN 
phase occurring in the late stages when the starburst and/or AGN essentially explosively 
drives out all the remaining interstellar material in the system. To test the predictions of 
the evolutionary sequence more robust SFR and AGN activity estimators are needed. This 
reinforces the importance of understanding the contributions of SB and AGN activity in 
known examples in which they coexist. 

Like most ULIRGs, IRAS 19254-7245, at z — 0.062, is a merger system of two gas- 
rich spiral galaxies, with the nuclei ~ 9" (which corresponds to ~11 kpc) apart from each 
other. Among these interacting ULIRGs, IRAS 19254-7245, with L Fm = 2 x 10 12 L Q , is 
a remarka ble one. It is kn own to host a powerful SB with star formation rate (SFR) ~ 
80M Q /yr (IBerta et al.ll2003l ) and a heavily obscured AG N, which resembles in many aspects 
the prototypical type-II (obscured) quasar NGC 6240 (jlwasawa fc Comastril Il998l ). This 
source exhibits giant tidal tails extending to a distance of about 350 kpc triggered by a 
merger event reminiscent of the cla ssical u Antennae" se en in NGC 4038/4039. The southern 



galaxy is classified Seyfert 2 (IMirabel et al. 



optical spectrum typical of SB (IBerta et al.ll2003l ) 



19911 ). and the northern galaxy shows an 



Previous X-ray observations of IRAS 19254-7245 are strongly indicati ve that this ULIRG 



harbors both a po werful SB and a deeply buried high- luminosity AGN (IPappa et al.l 12000 



Braito et al.l 120031 ) . A recent Suzaku observation has detected the intrinsic AGN emission 
in the 15-30 keV band and has confirmed the XMM- Newton detection of Fe-Ka emission 



(IBraito et al.ll2009l ). These data show that both neutral and ionized Fe-Ka lines are present. 
In this scenario, the intrinsic 2-10 keV luminosity of IRAS 19254-7245 could be greater than 
10 44 ergs s -1 , dominated by the buried quasar. Thi s picture has been confirm ed by high sig- 



nal t o noise L band (3 — 4/im) (IRisaliti et al.l 120031 ) and Spitzer observations ( INardini et al. 



20081 ) which also suggested that the main energy source in the mid-IR is an obscured AGN. 
The XMM-Newton and Suzaku observations have also shown that the soft emission is dom- 
inated by the starburst. The temperature kT ~ 0.85 keV and the luminosity of this soft 
component are consistent with the hot gas found in starburst galaxies and in the other 
ULIRGs observed so far. The XMM-Newton data show that, although the X-ray emission 
comes mostly from the southern nucleus containing the buried AGN, the source slightly 
extended. However, the XMM PSF was too broad to allow the structure of the soft X-ray 
emission to be investigated in more detail. Furthermore, since the smallest possible extrac- 
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tion region with XMM-Newton includes the whole system, it was not possible to spatially 
disentangle the contribution of diffuse SB and AGN emission. Again despite the intrinsic 
high luminosity of the AGN in IRAS 19254-7245, it is heavily obscured and thus it may not 
outshine the starburst emission below 2 keV. 

Here we present the analysis of three exposures of IRAS 19254-7245 by Chandra with 
total exposure time of 75 ks. These are the first X-ray observations that allow us to directly 
disentangle the AGN and starburst. We present the Chandra data reduction in Section 121 
and we analyze and fit the spectra of the different parts of IRAS 19254-7245 in Section El 
We then discuss our results and draw conclusions in Section HI 

2. Observations and data reduction 

The Chandra X-ray Observatory observed IRAS 19254-7245 for three separated expo- 
sures with the back-illuminated chip S3 of the Advanced CCD Imaging Spectrometer (ACIS) 
instrument in a FAINT mode in December 2009 (see Tabled]). The data were reprocessed 
using the contributed script chandra_repro of the CIAO software (ver. 4.3) and calibration 
database (CALDB 4.4.2) to create the new Level 2 event files. Sources were detected using 
the wavdetect tool. 

In Figure[U we show four images of IRAS 19245-7245 by merging the three observations, 
which are within the energy bands of a) 0.3-2 keV; b) 2-4 keV; c) 4-8 keV and d) 6-7 keV. 
The images are smoothed (in ds9) with a Gaussian kernel of 3 pixels (1.5" in radius, where 
1" corresponds to 1.2 kpc at the redshift of IRAS 19245-7245). We marked the southern (S), 
northern (N) nucleus and an off-nuclear point source (XI) in panel (b). As shown in panel 
a) in Figure [U there is extended soft X-ray emission contributed by the circumnuclear star 
formation regions, and we divided the diffuse emission into two regions in shape of annuli 
with their origin at the southern nucleus: Dl (3-10", i.e., 3.5-11.8 kpc) and D2 (10-20", i.e., 
11.8-23.6 kpc). In the following subsections, we will analyze the imaging and spectra of the 
nuclei and the diffuse emission. 

Throughout the work in this paper, all Chandra spectra as well as the associated re- 
sponse (RMF) and ancillary (ARF) files were extracted by running the script specextract, 
which has slightly different input parameters for point-like source and extended emission. 
We adopt cosmological parameters of H = 70 km s _1 Mpc -1 , Q M = 0.27 and fl\ = 0.73. 
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3. Analysis 



3.1. The entire merging system 



IRAS 19254-7245 was obs erved by XMM-N ewton in March 2001 (IBraito et all 120031 ) 
and Suzaku in November 2006 ( IBraito et al.ll2009f ). The published X-ray fluxes and spectral 
properties were derived from spectral analysis of the global emission (here 'global' means 
the entire region including the two nuclei and the extended starburst region). We thus first 
analyzed the global spectra of our Chandra observation, and compared them with the XMM 
and Suzaku observations to see if there is any variability. The radius of the extraction region 
is 15", centered on the southern nucleus, which follows the work of the XMM observation. 
The source region in the Suzaku paper has a radius of ~ 2.9', which covers nearly the entire 
S3 chip. There are several visually-detected point sources which are located 1-3' away from 
the southern nucleus of IRAS 19254-7245. Therefore, we did not use this source radius to 
extract the global spectra. The background emission was chosen from a close-by source-free 
circular region. 

The effective exposure time after flare filtering is shown in Table [TJ Photons with energy 
greater than 8 keV were ignored in the spectral fitting, since the effective area of the ACIS 
chips is very small above 8 keV. We obtained 1850 photons in total for all three exposures. 
The spectrum of each observation was grouped to at least 10 counts per bin to allow the use 
of x 2 statistics. 



Following the work of IBraito et al.l (120031 . 120091 ). we fitted each spectrum using a model 
of thermal plus power-law with an additiona l Fe Ka line. The weighted average Galactic H I 



column density is iVH,G = 5.95 x 10 cm 2 (jDickey &: Lockmanlll990l ; iKalberla et al.ll2005l ) 



and is set to be the lower boundary of the obscuration during the fit. We fixed abundance 
in the thermal component to solar. The photon index of the power-law continuum was 
determined by ignoring photons with energy greater than 6 keV and with no Gaussian 
component added in the model. We then fixed the photon index and used the photons in 
energy band 6-8 keV to fit the iron emission line. 

We plot the 0.5-2 keV and 2-10 keV fluxes of the Chandra and the other two observa- 
tions by XM M and Suzaku in Figure [2J The observed X-ray luminosities reported in XMM 
observation (IBraito et al.l 120031 ) were about a factor of two higher because of the cosmo log- 
ical parameters used in their paper, so we re-derived the luminosities by using the latest 
parameter values above. 

Since there is no significant variability of either the soft or hard X-ray flux between each 
exposure of the Chandra observation, it is possible for us to fit their combined spectra. The 
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spectra were summed using addspec in HEASof t FTOOLS. Again, the combined spectra were 
grouped for at least 10 counts per bin, and x 2 statistics was used in fitting. The best-fitting 
parameters are shown in Table [2J as well as the observed soft and hard X-ray luminosities 
contributed by each component, where errors are determined at the 90% confidence level for 
one parameter of interest. 

In order to put more constraints on the spectral properties of IRAS 19254-7245, we add 
the XMM and Suzaku data to the spectral fitting. We loaded the PN spectrum of XMM 
(as the first group), the XIS-FI and XIS-BI spectra of Suzaku (as the second group), and 
the three ACIS-S spectra of Chandra (as the third group) to perform a spectral fit of the 
time-averaged spectra. The parameters of the thermal and power-law components are tied 
together in different spectrum groups except for the relative normalization via a constant 
multiplicative factor. The temperature and power-law index parameters are close to the 
results of Chandra data only but with smaller uncertainties, as shown in Table [2j The 
luminosities are calculated as the mean of the three observations. The plasma abundance 
was initially fixed at solar, and we then allowed it to vary as a free parameter in this fit. 
The best-fitting abundance was Z = O.34Z , and was poorly constrained. We thus fix this 
parameter to be at solar in this and the following spectral analysis, since the model fitting 
is not very sensitive to it. We first tied the normalization of the iron line component in all 
three observations, and the line energy was found to be at 6.68 keV, consistent with the 
Chandra only fit above. However, the line EW was only ~ 0.78 keV compared to 1.71 keV 
we obtained above, which is a hint for variability of the iron line emission. Thus, we allowed 
the normalization of the iron line component to be free between the three groups. As the 
EW and luminosity of iron line differ between three observations, we leave them as blank in 
Table |2j and we will discuss the variability in the following sections. 



3.2. Southern nucleus 



The southern nucleus i s known to host a Compton-thick AGN by previous X-ray obser- 
vations (IPappa et al.ll2000l ; iBraito et al.ll2003l . 120091 ) . Those X-ray satellites have low spatial 
resolution, and the detected X-ray emission from the AGN was contaminated by the host 
galaxy. We are now able to disentangle the AGN from its nuclear starburst region in our 
Chandra observation. We extracted the spectra of the AGN in the southern nucleus from 
a circular cell with radius 2.5". The nucleus is a point-like source, so an energy-dependent 
point-source aperture correction was applied to the source ARF file. A local background 
spectrum was extracted from an annulus centered on the nucleus with an inner radius of 3" 
and an outer radius of 4". 
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We did a simultaneous joint fit to all three spectra with a model of an absorbed power- 
law plus a Gaussian at first. Then, an additional thermal component (APEC) was added to the 
fit. The abundance was again fixed at the solar value, and we found a best-fit temperature 
of the plasma kT = 0.79+°;| keV. x 2 decreases from 90.0 for 105 d.o.f (without the thermal 
component) to 81.7 for 103 d.o.f. Compared to the power-law plus Gaussian model, this 
model was preferred at 99.3% confidence level according to the F-test. 

We also fit the summed spectrum and the spectral parameters are close to the joint fit 
above. We list the best-fitting parameters of both the joint fit and the summed-spectrum 
fit in Table [2j The plot of the summed spectrum is shown in Figure [31 The flat continuum 
(r ~ 1.2) of the hard X-ray spectrum and strong iron emission line at 6.7 keV with EW ~ 
1.5 keV are consistent with the previous X-ray observations. 

An interesting question is if the soft flux in this spectrum is due to circumnuclear 
starburst emission or is more directly associated with the AGN. To test how this emission 
scales with extraction region size we also extracted the spectrum of southern nucleus in a 
smaller region with radius of 1.5". However, the background was extracted from the larger 
source region (with radius = 2.5", the same region as the source extraction region above), 
scaled to the area of the 1.5" circle. The rationale is that any component that is extended 
over the 2.5" region with an ~ constant surface brightness would subtract out completely, 
while nuclear point-like emission by reduced only by the ratio of the areas of the two regions, 
i.e., 1.5 2 /2.5 2 ~ 1/3. The temperature of the thermal component is kT = 0.78^'^, nearly 
the same as the value above but with larger uncertainty. The power-law normalization 
is two thirds of that in the spectral fit of the larger extraction region case, and so is the 
iron line normalization, consistent with both of these components being point-like. The 
ratio for the APEC normalization holds the same scaling relationship as that of the power- 
law normalization in the two cases. This indicates that either there is also point like, i.e., 
dominated by a thermal source in the innermost part of the nucleus or the AGN emission is 
mixed with the ambient star-formation and we cannot resolve it with Chandra. If the latter 
is true, we would expect consistency with the diffuse emission from the extended starburst 
region, which we will analyze in section 13.41 below. 



3.3. Northern nucleus 

Although the northern nucleus is relatively weak in X-ray emission compared with the 
southern one and is believed not to harbor an AGN, it is still interesting to investigate its 
X-ray properties. We positioned the northern nucleus by visual examination and extracted 
the spectra within a circle of radii of 1.5" for Obs ID 10538 and 11797. It could not be 
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visually detected in the observation of Obs ID 11798, and the spectrum was extracted from 
the region used in Obs ID 10538, which has the brightest northern nucleus in three exposures. 
30 photons in total are collected in 0.3-8 keV range. We summed the spectra a nd fitted th e 
combined spectrum with an absorbed power-law model using Cash statistic ( ICashl Il979l ). 
Although Cash statistics is designed for unbinned data, we groupe d the spectra to a t least 
1 count per bin, which could improve the performance of the fit (ITeng et al.ll2005l ). The 
photon index is found to be Y = 1.2^^. The observed 2-10 keV luminosity derived from 
this model is L 2 _iokeV = ^-0 X 10 40 erg s -1 . The spectral parameters are listed in Tabled 



3.4. Diffuse emission 

3.4-1- Imaging 

Thanks to the high spatial resolution of Chandra the diffuse emission from the star 
formation regions could be imaged (see Figured]). In order to study the extended structure, 
we investigated the brightness distribution of IRAS 19254-7245. We first simulated an 
observation of a point source wit h the same spectra l properties as the merged event using 



the Chandra Ray Tracer (ChaRT; ICarter et al.ll2003l ) to trace the rays through the Chandra 



X-ray optics, and created an image of the PSF using MARX (version 4.5). The comparison 
between the observed radial profile of the extended structure and the simulated PSF in both 
soft and hard X-ray bands is shown in Figure HI where the X-ray centroid of the southern 
nucleus was used as the origin. In the left panel of Figure HI we see that the observed 
distribution of the soft X-ray emission is rather close to the simulated one within the radius 
of ~2", and for radii between 3"and 20", there is a significant excess of the observed counts 
with respect to the simulated PSF. In the hard X-ray band shown in the right panel of Figure 
IU the radial profile of the extended emission is consistent with the simulated PSF. 



3-4-2. Spectral analysis 

In order to extract the spectrum of the circumnuclear starburst, the emission from 
the AGN was excised within a circle of 3" radius, and other point sources detected by 
wavdetect were also excluded. We extracted the spectra of the star- forming region in two 
annuli centered on the southern nucleus: one has the inner radius of 3" and outer radius of 
10" (Dl); the other one has inner radius of 10" and outer radius of 20" (D2), as described in 
Section [2j The background spectra were extracted in a concentric annulus with inner radius 
of 40" and outer radius of 50". The spectra of three exposures were then added to improve 
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the signal-to-noise ratio. 



Soft emission: The extended star-forming region is dominated by soft emission, as 
shown by the comparison between the soft and hard X-ray images in Figure [TJ We first 
ignored the hard X-ray photons with energy greater than 3 keV, and fitted the spectrum of 
the soft emission. It was initially characterized by a thermal plasma with solar abundance. 
The resulting temperature was kT = 0.77l°;°9 keV for Dl and kT = 0.91^ keV for D2 
(see Table [3]). The plasma temperatures are consistent with the spectral fits of XMM and 



Suzaku data, and a re within the the typical range of starburst galaxies (e.g. JPtak et al.lll999 



Grimes et al.ll2005l ). The temperature of the diffuse emission in Dl is also close to the value 
for the AGN spectra, indicating that the emission from the southern nucleus is probably 
also exhibiting thermal emission from the circumnuclear star-forming regions. We then 
investigated the possibility of introducing another thermal component assoc iated with an 



absorber in the spectral fitting since such a feature is fo und in other ULIRGs (IDahlem et al. 



19981 ; llwasawa fc Comastrilll998l ; iGallagher et al.ll2002l ). This second plasma component is 
colder, with a temperature of kT = 0.34^ keV for Dl and kT = 0.27±%{l keV for D2. 
The best-fitting parameters are listed in Table [3J The two-temperature model does not 
improve the spectral fit much, and the normalization of the cooler component has a very 
large uncertainty, so we thus adopt the results of the single temperature model in this paper. 

Hard emission: The hard continuum was fitted by a power-law, and we included the soft 
emission in the spectral fitting because the power-law parameters could be determined much 
more accurately in this way. For the inner region Dl, the photon index is 1.73±°J°, a typical 
value for AGN. For the outer region D2, we initially got a much steeper index of 2.70lg §g. 
However, the larger errors due to low signal-to-noise makes it consistent with the value in 
the spectral fit of Dl. The soft and hard X-ray luminosities are also listed in Table [2], and 
the spectral plots of both regions are shown in Figure As we see in Table [21 the thermal 
emission and power-law emission account for similar amounts of the soft X-ray emission. The 
hard X-ray luminosity is 20 times smaller than that of the southern nucleus. The best-fit 
column density along the light of sight is negligible in D2 so we set it to be equal to the 
Galactic column density, and it is slightly larger in Dl. Because the hard X-ray emission in 
the D2 is much weaker than in Dl and has poorer quality, we extracted the spectra from the 
combined regions of Dl and D2 and re-fitted the hard spectrum. The power-law index was 
then found to be close to that in Dl. 

The spectrum of the diffuse emission shows an excess in the 6-7 keV band. Could 
i t be due to ionized iron emission contributed by very hot gas driven by the starburst 
(IHeckman et al.lll990l ; IStrickland fc Heckmanl 120071 )? We show in panel d) in Figure [H the 
image within the energy range of 6.0-7.0 keV, which covers the neutral and ionized iron 
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emission line energies in observer's frame. This image shows that the region of iron line 
emission overlaps the southern nucleus and is not extended. However, due to the low signal- 
to-noise ratio and limited photons in hard X-ray band, we cannot rule out the possibility 
that the hard continuum emission is due to another much hotter plasma. If we fit the hard 
X-ray spectrum with a thermal component, the temperature is kT ~ 6 keV with very large 
uncertainty, and the reduced x 2 is nearly the same as that in the case of power-law. 



3.5. Iron line emission 

The i ron emission line in the Southern nucleus has a relatively broad profile (a ~ 0.3 



keV), and lBraito et al.l (120091 ) showed the possibility that it is a blend of multiple lines. It is 
therefore interesting to further investigate this feature. In order not to wash out the possible 
fine structure of the iron line profile, we fitted the unbinned time-averaged Chandra spectra 
using the Cash statistic. 

We used a power-law continuum with the photon index fixed at T = 1.22 and two 
gaussian lines to fit the Chandra spectrum within energy between 3 keV and 8 keV. The 
centroid energy of two lines are found to be 6.43 and 6.71 keV, and the equivalent widths 
are 0.17 keV and 0.66 keV respectively, i.e., the neutral iron line is weaker than the ionized 
iron line. However, this is not the same as when we re-analyzed the XMM data. We fitted 
the unbinned PN data of the XMM observation in the same way as above. The energy of 
the ionized iron line is 6.89 keV, which is Hydrogen-like. The EW is 1.16 keV for the neutral 
line, and 0.44 keV for the ionized line. The double line fit of Suzaku data gives similar results 
as Chandra for the line energy and EW as the Chandra observation. We then plotted the 
confidence contours between the line energy and normalization for both lines at the la, 2a 
and 3<T levels. The line profile and contours are shown in Figure |6l where the energy is 
in the rest frame. We see that the significance of the 6.4 keV line in Suzaku and Chandra 
observations is only 2cr, but the lines in the XMM observation could be separately detected 
at > 99% confidence. 

We then ran Markov Chain Monte Carlo (MCMC) simulations to test the statistical 
significance of the lines. Each chain entailed 10 5 iterations, and we ran 5 chains for each 
spectrum. The parameters and the associated errors are given in Table HJ The significance 
of the line was then calculated as the fraction of the iterations where the line normalization 
parameter was non-zero. We see in Table H] that the 6.4 keV line in the XMM observation 
and 6.7 keV lines in the Suzaku and Chandra observations are significant at a level of over 
99.99%. The significance for the 6.9 keV line in the XMM observation is about 97%, and 
only less than 90% for the 6.4 keV lines in the Suzaku and Chandra observations. As in 
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Figure [6], we plot the probability contours of line energy and line normalization derived from 
the MCMC analysis in the upper panel of Figure [7] at probability levels of 4 x 10~ 4 , 8 x 10~ 4 
and 1.2 x 10~ 3 . Histograms showing the distribution of each line normalization is also shown 
in the lower panel of Figure [71 The left plot shows the neutral 6.4 keV line normalization 
distribution and the right plot shows the ionized line (6.7 or 6.9 keV) distributions. Lines 
in black, blue and red indicate the observation of XMM, Suzaku and Chandra, respectively. 
The peak of distribution of the 6.4 keV line overlaps in the Suzaku and Chandra observation, 
but is obviously separated from the XMM observation, giving a ~ 90% probability of the 
line variability. However, the variability of the ionized iron line is not significant as indicated 
by the right panel of Figure [7J 

We also tried to use three lines to fit the iron emission spectra. However, the data 
quality does not allow us to investigate the possibility that both He-like and H-like lines are 
present. 



3.6. Off-nuclear point source 



Ultraluminous X-ray objects (ULXs) are commonly seen in ULIRGs (IPtak et al.ll2003l ). 
They are objects that are much more luminous than a stellar mass black hole or neutron 
star binary accreting at Eddington limit, and are not located in the nucleus of a galaxy. As 
shown in Figured], an off-nuclear point source (XI) is found about 8" south of the southern 
nucleus, which is also detected by wavdetect. No counterpart at this position is found in 
optical images. We plot the light curve of XI in Figure |8] within the three Chandra exposures, 
and the bin size is 2 ks. The spectra of the non-nuclear point source were extracted from a 
circle of 2" radius, and the background region was chosen to be a concentric annulus with 
an inner radius of 4" and an outer radius of 5". This point source has 78 photon counts in 
the 0.3-8 keV band detected in three exposures. Thus we binned the summed spectrum into 
at least 1 count per bin using Cash statistics to perform the spectral fit. 

We fitted the spectrum of this point source with several models: simple power- law (PL), 
a thermal plasma (APEC) plus a power-law and a multi-color accretion disk (MCD) with 
a power-law. Table [5] lists the best-fitting parameters and the observed soft and hard X- 
ray luminosity of each model. The best-fitting photon index is about T = 2 in the PL and 
APEC+PL model, and a little bit harder (r = 1.6) in the MCD+PL model. The temperature 
of the plasma in the APEC model is found to be kT = 0.3 keV, and the temperature at 
inner disk radius of the MCD model is kT- m = 0.46 keV. All three models give the observed 
soft X-ray luminosity of about 2.4 x 10 40 ergs s _1 , and the hard X-ray luminosity of about 
4 x 10 40 ergs s _1 , which is about one order of magnitude more luminous than a typical ULX. 
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4. Discussion 

The Chandra observation allows us to disentangle the two nuclei in the merger system 
of IRAS 19245-7245 for the first time in X-ray imaging. Our analysis demonstrates that 
the southern nucleus, which hosts a Compton-thick AGN, contributes over 80% of the X-ray 
emission of the entire merger system in 0.5-10 keV band, and produces ~ 90% of the hard 
X-ray luminosity. Extended soft X-ray emission has been observed around the southern 
nucleus out to a radius of ~ 25 kpc, contributing a similar amount of soft X-ray photons as 
the southern nucleus itself. 



4.1. Emission from AGN 



Soft emission: The flat hard X-ray contin u um and th e intrinsic Compton-thick ab- 
sorber have been well discussed by iBraito et al.l ( 120031 . 120091 ) . One of the key goals of our 
Chandra observation was to better investigate the relative contribution of SB and AGN 
activity. However, this discrimination becomes vague when referring to the soft emission 
from the southern nucleus. The soft emission is fit well with a power-law and a thermal 
component. The properties of the thermal component are consistent with those from the 
extended region and thus is likely to be thermal and starburst dominated emission rather 
than AGN-photoionized plasma. Even when we extracted the emission from a smaller region, 
this thermal component could not be ruled out. The total observed 0.5-2 keV luminosity of 
the southern nucleus is 2.6 x 10 41 ergs s _1 , and one third comes from the thermal emission 
according to the spectral model. If we conve rt the soft X-ray luminosity to SFR assuming 
the relationship given by iRanalli et al.l ( 120031 ). we get SFR = 57 M Q yr" 1 , which is close to 
the estimated SFR that is about 80 M Q yr -1 for the southern nucleus ( jBerta et al.l 120031 ) 
given the dispersion in the X-ray/SFR correlation. 

Fe Ka line: Iron emission lines are only detected in the southern nucleus. We tried 
to add a Gaussian component in the spectral fitting of the inner region Dl in the extended 
starburst region, however, it did not reduce the x 2 - The low quality of the hard X-ray photons 
in the outer region D2 prevents us from investigating the iron line feature. Indeed, the 6-7 
keV image shows that iron line photons are only from the southern nucleus. The luminosity 
of the iron emission line observed by XMM-Newton in March 2001 is 3.3 x 10 ergs s _1 (after 
K— correction), and it dropped to 2.2 x 10 41 ergs s _1 in the Suzaku observation in November 
2005. The result from our recent Chandra observation is, however, consistent with the XMM 
observation. The EW of line was 1.5 keV and 0.7 keV in XMM and Suzaku observations, 
respectively. 
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We also investigated the possibility of a blend of two lines of the Fe Ka emission. 
In the 2001 XMM observation, a neutral iron line with centroid line energy at 6.4 keV 
and an ionized H-like iron line with centroid line energy at 6.9 keV are clearly separated 
at significance level > 99%. The intensity of the 6.4 keV line is larger than the 6.9 keV 
line. However, the significance of the 6.4 keV line in the 2005 Suzaku and 2009 Chandra 
observations decreases to only 2a, and a He- li ke iron line (6.7 keV) sh ows up in the late 



two observations and dominates the intensity. iTerashima et al.l (120001) reported a s imilar 
iron line variability in NGC 4579 in two ASCA observations, but lEracleous et al.l (120021 ) 
compared that with a Chandra observation and argued that such variability might be due to 
the very large extraction region of ASCA spectra, which was contributed by the ionized Fe 
Ka emission from the entire circumnuclear nebula. However, as we show in above sections 
that there is little extended emission in 6-7 keV band, the aperture effect is not likely the 
cause of the variability. 

The variability of the 6.4 keV line indicates that the cold matter where this line originates 
lies within 2 pc from the nucleus. If the energy shift of the ionized line from 6.9 keV to 6.7 
keV is true, there might be multiple locations of the clumpy clouds with different ionization 
states. Since there is no flux variability of the continuum emission, we could not constrain 
the distance of the ionized line emitters based on the time lag between the continuum and 
the Fe line ionization change. If the change of ionization state (from H-like to He-like) is not 
true but due to the insufficient photons in the spectral fitting, there is still a ~ 70% change 
in the ionized Fe Ka line intensity. The 6.7 keV line could originate either from an ionized 
accretion disk or from reflection/scattering of the hard nuclear X-rays by ionized clouds not 
farther than 2 pc from the central engine. 



4.2. Emission from Starburst 

Soft emission: From the soft X-ray image of IRAS 19245-7245, we can see that hot gas 
emission extends to a radius of ~ 20" (~ 25 kpc) surrounding the southern nucleus. The 
thermal emission requires a kT ~ 0.8 keV component in the best-fitting models of the inner 
region (Dl) of the diffuse emission. The temperature is consistent with the thermal emission 
found in the southern nucleus. Although the spectra in the outer region (D2) have relatively 
poorer quality, the temperature of the hot gas in D2 from the spectral fit is consistent with 
that of Dl but with larger uncertainty. 

Using the parameters of the X-ray spectral model, we can derive the basic physical 
properties of the hot gas in the regions of S, Dl and D2. The normalization for the APEC 
component (10 -14 J n e nndV / 'Atc\D ^(1 + z)] 2 ) implies a volume integral of density squared if 
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we assume the electron density (n e ) approximately equals the H atom density (riff), where 
D A = 243.3 Mpc is the angular distance of IRAS 19254-7245 at redshift z = 0.062. The 
mass (M), mean thermal pressure (P), total thermal energy (E) of the hot gas can be 
derived using the electron density and temperature from the spectral fit. The characteristic 
timescale could be calculated as the size divided by the sound speed in the hot gas. Then, 
the outflow rates of the mass (M) and energy (E) can be roughly estimated according to 
the timescale. These physical parameters are listed in Table El where we parameterize them 
by the volume filling factor / of the X-ray-emitting material. We see that while most of the 
mass and thermal energy content of the hot gas is contained in the outer regions (Dl and 
D2), the actual estimates of the rates at which energy and mass flow outwards are consistent 
as we move from the central region (S) to Dl and through D2. This is constituent with a 
mass- and energy- conserving outflow which originates in the starburst (inside region S). 



We then use Starburst99 (ILeitherer et all Il999l ) to estimate the relevant parameters 
expected for the collection of sup ernovae and massi ve stars in the starburst. We adopt a 
star formation rate of 80 M Q yr -1 ( iBerta et al.ll2003l ) which has been constant over the past 
100 Myr and has a standard Kroupa IMF. The resulting total mechanical power of the stellar 
wind and supernovae ejecta is 4 x 10 43 erg s _1 , they carry a momentum flux of 3 x 10 35 dyn, 
and the mass outflow rate is Msn ~ 16 M Q yr _1 . 

These values imply that a starburst- driven outflow can account for the properties of the 
extended soft X-ray emission. The observed outflow rate of (thermal) energy is only 10% 
of the total energy injection rate from the starburst while the observed mass outflow rate is 
about three times larger than the expected starburst mass injection rate. This latter result is 
commonly observed in starburst-driven winds and it ascribed to ma ss-loading as the super- 
novae and massive star ejecta mix with ambient interstellar gas (e.g. IStrickland fc Heckman 
20091 ). The predicted momentum flux from the starburst results in a wind ram pressure of 
about 1.8 x 10~ 10 , 1.6 x 10~ n , and 4 x 10~ 12 dyne cm -2 at the outer boundaries of regions S, 
Dl, and D2 respectively. These values are very close to the inferred pressures in the observed 
extended soft X-ray region. 

More empirically, we can compare the basic properties of the soft extended emission 
to what is seen in typical starbursts. First, we can compare the relation between the far- 
IR and soft X-ray luminosity of IRAS 19245-7245 with other star-forming galaxies which 
differ in their size and luminosity. The plot is shown in the upper panel of Figure [9j where 
the other galaxies are dwarf starbursts, n ormal starbursts, Lyman b r eak analogs, starburst - 
dominated ULIRGs and AGN ULIRGs flGrimes et all I2OO5I . bood . I2OO7I : Ijia et al.l lioilh . 
The soft X-ray luminosities are derived from the thermal component in the spectral fits by 
excluding the contribution from AGNs and X-ray binaries. The luminosity ratio of SX/IR 
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for ULIRGs is found to be about one dex smaller than the value for star-forming galaxies 
with Lir < 1O 11 L . This deficiency in systems with h igh IR luminosity was c haracterized 



as < logL S x/£iR >= —4.55 with la scatter of 0.37 by ISymeonidis et al.l ( 1201 ll ). This ratio 
for IRAS 19245-7245 is -4.66 according to our spectral analysis, well consistent with their 
result for the star-forming galaxies with similarly high far-IR luminosities. We also show in 
the lower panel of Figure |9]that IRAS 19245-7245 obeys the relation between the size of the 
soft X-ray emitting region and the far-IR luminosity defined by other star-forming galaxies. 

We conclude that the properties of the large (~ 50 kpc diameter) spatially-extended 
region of soft X-ray emission are entirely consistent with an origin as a starburst-powered 
galactic wind. 

Hard emission: The hard X-ray continuum of the diffuse emission is characterized as a 
power-law with photon index of T ~ 1.8, which is much steeper than that found in southern 
nucleus. We compare the luminosity ratio of hard X-ray to far IR with other LIRGs and 
ULI RGs observed by Cha ndra in F igure [TOl The o ther samples are from llwasawa et al.l (120091 ) 
and iLehmer et al.l (120101 ). where llwasawa et all s sample was divided into two groups of 
AGNs and "Hard X-ray Quiet" galaxies (HXQs) indicated by diamond and asterisk symbols, 
respectively, and ILehmer et all s sample was indicated by triangle symbols. We show IRAS 
19254-7245 as plug signs in two positions in this plot: one with AGN contribution included 
(in red) and the other one without AGN contribution (in blue). We see that the hard X-ray 
luminosity produced by the star formation process is consistent with other pure star-forming 
galaxies (with no AGN) having similar far-IR luminosities. 



4.3. Emission from other regions 



Northern nucleus: The northern nucleus was not resolved in the XMM-Newton observa- 



tion, and it was included in the 15" extraction circle of spectral analysis (IBraito et al.ll2003l ). 
According to our Chandra observation, both the soft and hard X-ray luminosities are ~ 30 
times smaller than those of the southern nucleus. We fitted the spectrum of northern nucleus 
with only an absorbed power-law due to the limited photons we obtained. The photon index 
is found to be T = 1.25 but poorly constrained. It has been demonstrated by the optical 
spectra that there is no AGN emission in this nucleus. This classification is supported by 
its low X-ray luminosity and the non-detection of Fe Ka line emission in our Chandra ob- 
servation. The SFR of the northern nucl eus is only ~ 0.5 M yr -1 , and it appears to be in 
a post-starburst phase (IBerta et al.ll2003l ). which is consistent with its weak X-ray emission. 



Off-nuclear point source: One off-nuclear point source is found ~ 8" south of the south- 
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ern nucleus. If we assume that it is associated with the Superantennae, it would lie about 
9.4 kpc from the southern nucleus in projection and have a 0.5-10 keV luminosity as high as 
6 x 10 40 ergs s _1 (based on the spectral fits using different models). This is about one order 
of magnitude more luminous than the typical ULXs found in other galaxies, but still within 
the lumi nosity rang e they span. In a recent Chandra survey of X-ray point sources in nearby 



galaxies, [Liu] (120111 ) found 479 ULXs with L x (0.3 - 8 keV) > 2 x 10 39 erg s" 



and about 

15% (74 out of 479) are extreme ULXs with L x (0.3 - 8 keV) > 10 40 erg s -1 . The X-ray 
luminosity of this ULX in the Superantennae would imply a minimum black hole mass of 
about 500 M© for a source at the Eddington limit. 



We note that an association with the Supernantennae is likely. Using the 2-8 keV flux 
of ~ 4 x 10~ 15 erg cm" 2 s -1 and the numb er-count relations from Chandra Deep Fields 



(IBauer et al.l I2004J ; iBrandt fc Hasingerl 120051 ) . we calculate that the possibility of a back- 
ground AGN located within a radius of 8" from the southern nucleus is only about ~ 1. 



o. 



5. Summary 

In this work we have presented Chandra observation of the proto- typical merger/ULIRG 
IRAS 19254-7245 (the Superantennae). This is the first X-ray observation that could spa- 
tially resolve the diffuse emission of this ULIRG. The main results are summarized as follows: 

1. The southern nucleus hosts a Compton-thick AGN, which contributes ~ 90% of the 
hard X-ray emission of IRAS 19254-7245. The soft X-ray emission from the AGN 
region is contaminated by the circumnuclear starburst. 

2. The diffuse emission from starburst extends out to a radius of about 25 kpc surrounding 
the southern nucleus, and is dominated by thermal emission with a temperature of 
~ 0.8 keV. This emission accounts for ~ 33% of the soft emission from the whole 
merging system. 

3. The global properties of the extended thermal emission are consistent with winds seen 
in many starbursts and the basic physical properties are fully consistent with an outflow 
driven by the supernovae and massive stars in the starburst. Feedback from the AGN 
is not necessary to explain what we see. 

4. The prominent Fe Ka emission line is a blend of both neutral and ionized lines, and 
variability is detected for the neutral one. 

5. No evidence for an AGN in the northern nucleus is detected in this Chandra observa- 
tion. 
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6. An off- nuclear point source is found and could be classified as a ULX if it is associated 
with the Superantennae. 
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Fig. 1. — The Chandra images of IRAS 19245-7245 in the energy bands of a) 0.3-2 keV; 
b) 2-4 keV; c) 4-8 keV; and d) 6-7 keV. North is up, and east is left. These images are 
merged from the three separated exposures, and smoothed using a Gaussian function with 
kernel radius of 3 pixels. The annuli of Dl (3-10") and D2 (10-20") are indicated by the green 
circles, and the das hed white circle indicates the extraction region (15") in XMM observation 
(IBraito et al.l 120031 ) . The southern (S) and northern (N) nuclei are marked in panel b), as 
well as one off-nuclear point source (XI). 
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Fig. 2. — The soft (asterisk) and hard X-ray (plus) flux of each observation. The points cor- 
respond in chronological order to XMM, Suzaku and Chandra observations. X-axis indicates 
the days since the date of XMM observation, and is not in linear scale. 
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southern nucleus 




Energy (keV) 

Fig. 3. — Chandra spectrum of the southern nucleus. The spectral model consists of a 
thermal, a power- law and an emission line component. Photon energy is de-redshifted to the 
rest frame. Top, data and the model; bottom, ratio of data to model. 
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Fig. 4. — Radial count distributions of the observed extended structure from the centroid of 
the southern nucleus (in black) and the simulated PSF (in red) in both soft (left) and hard 
(right) X-ray bands. 




Fig. 5. — Spectral fits of the diffuse emission in annuli of 3-10" (left) and 10-20" (right), 
where the southern nucleus is at the center. 
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Fig. 6. — Upper. Iron emission line profile fitted by double Gaussian lines centered for the 
XMM, Suzaku and Chandra observations. The energy is de-redshifted to the rest frame; 
Lower. Confidence contour plots for the double Gaussian lines profile. Contours correspond 
to the 68%, 90% and 99% confidence levels respectively. 




Fig. 7. — Upper. Probability contours of the iron lines from MCMC. The contour levels 
correspond to probabilities of 4 x 10 -4 , 8 x 10~ 4 and 1.2 x 10~ 3 ; Lower. Distribution of 
the normalization parameter of the iron line component at 6.4 keV (left) and 6.7 or 6.9 keV 
(right) in XMM (black), Suzaku (blue) and Chandra (red) observations. The parameter 
probability is in arbitrary units, and the relative scale depends on the bin size. 
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Fig. 8. — The light curve of XI with bin size of 2 ks obtained from the three Chandra 
exposures. The time axis is broken for display purpose. 
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Fig. 9. — The total thermal X-ray luminosity (upper) and the radius enclosing 90% of the 
total flux in the 0.3-2.0 keV X-ray band (lower) vs. the far-IR luminosity. IRAS 19254-7245 
is marked in red cross symbol. Other sources in this plot are dwarf starbursts (square), 



starbursts (asterisk), Lyman break analogs (plus ), starburst-domina t ed ULIRGs 
and A GN ULIRGs (circle), which are taken from lGrimes et all J2005i bood . l2007h 
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Fig. 10. — Hard X-ray luminosity vs. far infrared luminosity. IRAS 19254-7245 is marked 
by plus sign in red (entire galaxies in cluding th e AGN) and in blue (without the AGN). 
Others are from llwasawa et al.l (120091 . Iw09) and iLehmer et al.l (12010L LehmerlO) for com- 
parison, where the diamond symbols indicate the Iw09 AGN samples, and asterisks and 
upper limit arrows are "Hard X-ray Quiet" galaxies (HXQs) from Iw09, and LIRG samples 
from LehmerlO are marked as triangles 
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Table 1: Observation log. 



Obs ID Date 


Exp. Time a 


YYYY-MM-DD 


(ks) 


10538 2009-12-14 


29.06 


11797 2009-12-19 


35.16 


11798 2009-12-19 


11.81 



a The exposure time is the good time interval after flare filtering. 



Table 2. X-ray spectral model parameters. 



Region 7V H a fcT b T c E Fc d a c EW f x 2 /" £sx,th g ^sx, P i h W W 



Global k 


1 


04+2.75 


0. 


79+0.12 
■ ' z -0.14 


1 


on+0.18 
ou -0.16 




0.13: 


-0.17 
-0.09 


1 71+0.90 
1 -' i -0.74 


146.0/154 


1 3+ 2 1 


9 7+0.8 

z -'-0.7 


22.6+™ 


4 1+ 2 - 1 

^-J— 1.8 


Global 1 


1 


oq+o.77 

zy -0.45 


0. 


77+O.O6 
•' '-0.07 


1 


n R +0.09 
zo -0.07 


6.68ir 5 5 


0.17: 


H0.06 
-0.05 




388.9/361 


1 7+O.8 
1 -'-0.4 


3 6 +0 - 4 

O.D_ 2 


22.7± 2 ; 4 




S (joint) 


3 


QC+1-57 
.oo_ 1 22 


0. 


7O+0.39 
• ' y -0.23 


1 


9O+0.25 
■ zo -0.17 


6-68±8:i? 


0.23: 


HO. 50 
-0.23 


1 - OZ -0.73 


81.7/103 


7 +0A 


9 1+O.6 
z - 1 -0.2 


19.4 + |f 


o + 1.4 
0.0_ 1 2 


S (summed) 


5 


9 R +2.87 
• zu -2.63 


0. 


7Q+0.29 
' y -0.19 


1 


99+0.19 
• zz -0.21 


6-69±8:?l 


0.24: 


HO. 29 
-0.13 


1 - ou -0.60 


99.1/97 


9+ 1 ' 1 

u - y -0.5 


1 7+O.3 
1 -' -0.3 


19.8+11 


q q+1.0 
°- J -1.3 


N 


2 


4Q+7.04 
■^ y -2.49 






1 


9 r+0.93 
— 1 23 










25.5/29 




10+ 005 

u - iu -0.04 


0.48lg;S 




Dl 





fiO+0.94 
• Oy -0.69 


0, 


78+0.13 
•'°-0.07 


1 


70+0^60 
•' J -0.56 










23.1/26 


U.OO_ 24 


u - + -0.18 


0.88l°i? 




D2 






0, 


OQ+0.46 
yu -0.27 


2 


70+0.87 
•' u -0.88 










28.5/28 


o-i7±8:ii 


25+° 15 

u - zo -0.17 


o.i5±8:S8 




D1+D2 


1 


■ 1U -1.02 


0. 


77+O.O8 
•' '-0.11 


1 


70+0-75 
'' y -0.59 










53.8/55 


o.86lS:£ 


n 4O+0- 35 

u -^ y -0.22 


1 1 c+0.82 
1 - lo -0.52 





Note. — Uncertainties correspond to the 90% confidence level for one interesting parameter. 

a Column density in units of 10 21 cm~ 2 . 

b Temperature of thermal plasma in keV. 

c Photon index of hard continuum. 

d Energy of line center in keV in rest-frame. 

°Emission line width in keV. 

Equivalent width of emission line in keV. 

s Soft X-ray luminosity contributed by thermal plasma in units of 10 41 ergs s _1 . 
h Soft X-ray luminosity contributed by power-law continuum in units of 10 41 ergs s _1 . 
'Hard X-ray luminosity in units of 10 41 ergs s _1 . 
J Iron emission line luminosity in units of 10 41 ergs s _1 . 
k Using stacked spectral data from Chandra observation. 
'Using Chandra, XMM and Suzaku data. 
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Table 3. Spectral fits parameters of the diffuse emission in 0.5-3 keV. 



Region 




kT x 


K* iV H , 2 kT 2 


K* 


x 2 /^ 




10 21 cm- 2 


keV 


10 21 cm- 2 keV 








Single temperature 


Dl 


N H ,G 


n 77+O.O8 
u - ' ' -0.09 


4.3+p; 2 x 10~ 6 






19.15/22 


D2 


^h,g 


fl Q1+°-48 
U - yi -0.60 


1.2±^ g x 10~ 6 






18.9/17 


Double temperatures 


Dl 
D2 


zi - y -10.2 

0.75+ 7 - 25 


u - J *-0.06 

n 27+ - 18 

U - Z '-0.12 


3.5l 3 y u x HT 4 TVh.g 0.76™ 
2.4™ x 10- 6 3.54±|;iS 0M±° Q f 4 


4 .8ir 

4-o±S:S 


x 10~ 6 
x 10~ 6 


15.2/19 
18.1/13 



*Plasma model normalization in units of 10 _14 /{4tt[Z)a(1 + 2)] 2 } J n e riYidV, where Da 
is the angular distance, and n e and nn are the electron and hydrogen number densities, 
respectively. 



Table 4. Iron emission line features. 











Spectral fit 












line 1 








line 2 






Observation 


E 




EW 




E 


a 


EW 




XMM 

Suzaku 

Chandra 


6.39™ 
6 4S+ ' 37 


o n+ - 27 

0.01 (fix) 
0.1 (fix) 


1 ifi+ - 94 

i - iD -0.50 

0.09™ 

0-17±8;i; 




6 89+ - 07 
6 69+ 04 

D - oy -0.04 

6 71 + - 07 

D -' i -0.03 


0.01 (fix) 

o.iol8:Sg 
o.o4i8:8S 


0-44™ 

0.4818:1? 

0.66™ 












MCMC 












line 


1 






line 


2 




Observation 


E 


a 


EW 


significance 


E 


(j 


EW 


significance 


XMM 

Suzaku 

Chandra 


fi 41 + - 10 
fi 41+° 18 
fi S5+ 21 


is+ - 07 

u - io -0.08 

0.01 (fix) 
0.01 (fix) 


i ir+0.54 
1 - lo -0.48 

0Q+ 008 

u - uy -0.09 

14+ 017 


> 99.99% 
87.79% 
81.16% 


6 8Q+ - 10 

°- oy -0.09 
r ro+0 06 

6 71 +0- 06 

D -' i -0.05 


0.01 (fix) 

o.ni8:S8 
o.o8l8:8 


n 44+°- 32 

U '**-0.34 

0.4718:1s 
0.89™ 


97.14% 

> 99.99% 

> 99.99% 



Note. — E, a and EW are in units of keV. 



Table 5. Spectral fits parameters of the off-nuclear point source (XI). 



model 


Nu 


kT or fcTin 


r 


C-statistic 


£o.5-2kcV 


i2-10kcV 




xlO 21 cm" 2 


keV 




/dof 


xlO 40 ergs s" 1 


xlO 40 ergs s" 1 


PL 

APEC+PL 
MCD+PL 


1 07+0. 81 
1 -°'-0.67 

2 7Q+0- 79 

z - ' y -0.76 
A+LG 


o.29l8:i 8 
o.46l8:S 


O pw+0.82 

z.uo_ 65 
2 07+ ' 77 

z - u ' -0.71 
1 - oo -0.30 


50.0/55 
49.0/53 
49.7/54 


9 4 +o.5 

9 0+O.6 
z -°-0.5 

2.4l?i 


4 - 1 -0.8 

4 0+ 1 ' 1 
4 5+ 1 ' 4 
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Table 6: Physical properties of the diffuse hot gas in regions S, Dl and D2. 



Parameters S Df D2 



n e (10 2 cm 3 ) 


5.5/ " 1/2 


0.44/~ 1/2 


0.087/ 2 ~ 1/2 


M (10 8 M ) 


2.5/ 1/2 




10.3/ 2 1/2 


M (M yr~ l ) 


39/ 1/2 


47/^ 


48/ 2 1/2 


E (10 56 erg) 


8.7/ 1/2 


25A 1/2 


36/ 2 1/2 


E (10 42 ergs" 1 ) 


4.3/ 1/2 


5.2/y 2 


5.3/ 2 1/2 


P (lCT 10 dyn cm" 2 ) 


1.6/ " 1/2 


0.13/^ 1/2 


0.025/ 2 _1/2 



Note. — fo, fi and fa are the volume filling factors of the X-ray emitting material in regions of S, Dl 
and D2, respectively. See text for details. 



